Specific Heat of the Periodic Anderson Model: Prom Weak to Strong Coupling 
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We study the temperature dependence of the specific heat in the periodic Anderson model as 
function of the on-site Coulomb interaction, hybridization, and position of the /—electrons energy 
level. At strong coupling (17 = oo) we use slave bosons, whereas at weak coupling we use a simple 
Hartree-Fock decomposition. We find that both in the strong and weak coupling limits the specific 
heat of the system presents a multiple peak structure in a mean field treatment. We believe these 
features to be related to band-shape effects. The temperature evolution of the low-temperature peak 
position on the model parameters for the non-half filled non-symmetric case is discussed. 

PACS numbers: 75.20.Hr, 71.27. +a, 71.28.+d, 65.40. +g 



Heavy-fermion materials are compounds containing 
rare earth and actinide atoms that exhibit high specific 
heats and spin susceptibilities at low temperature. This 
feature is associated with high effective masses of quasi- 
particles due to strong Coulomb interactions among the 
4/— or 5/— electrons.!!] The behavior of these materials is 
controlled by three parameters: the hybridization V be- 
tween /— and d— orbitals, which confers the /—electrons 
an itinerant character, the value of the /—level eo energy 
relatively to the chemical potencial, and the value of the 
Coulomb energy t/.El 

The strong Coulomb interaction together with the 
small values of the hybridization matrix allow for for- 
mation of local moments leading to magnetic behavior 
in these systems. The Kondo effect and the magnetic or- 
der due to the RKKY (Ruderman-Kittel-Kasuya-Yosida) 
interaction are consequences of the iaterplay between 
Coulomb interaction and hybridization.Eru Together with 
magnetism, heavy-fermion materials may present super- 
conducting ground states, whose .origin and nature are 
far from being fully understood.clEl 

The periodic Anderson model (PAM) is considered a 
good candidate for describing heavy-fermion systems. It 
is generally accepted that a complete understanding of 
heavy-fermion system properties requires going beyond 
mean field theory. Nevertheless, mean field theory still 
gives qualitatively correct information about some prop- 
erties of heavy-fermion systems J3 

Recently, the temperature dependence of the specific 
heat of the PAM, at half filling, was studied using the 
equation of motion method, both in the symmetric and 
non-symmetric cases,™ and a multiple peak structure 
was found. This multiple peak structure was first ob- 



tained by Zielihski and Matlak in the half filled symmet- 
ric (2eo + U = 0) case.£3 However, it is believed that 
mean field theory is not able to capture this multiple 
peak structure in the specific heat. 

This belief contrasts with Monte Carlo results show- 
ing that the magnetic Hartree-Fock approximation repro- 
duces well the free .energy of the PAM at low tempera- 
tures when V <C i7,Eil implying that the thermodynamics 
must be well described by the Hartree-Fock approxima- 
tion in this limit. 

In this study we report results for the PAM's spe- 
cific heat both in the weak and infinite coupling limits. 
Our study is done away from half filling and in the non- 
symmetric case. The strong-coupling limit U — oo is 
studied using slave bosons, whereas the weak coupling 
limit is treated within Hartree-Fock approximation. We 
find that in both limits the specific heat presents a multi- 
ple peak structure and that its dependence on the model 
parameters is in qualitative agreement with what is phys- 
ically expected. 

The standard periodic Anderson model model for spin 
1/2 electrons is written as 

H = H%+H° f + H df +Hu, (1) 

where 

H f = E( £ " - Mich' > ( 3 ) 

Hdf=VY J {4,Ji,<r + fiji,v), (4) 
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and eg is the dispersion of the d— electrons, eo is the bare 
energy of the localized /—states, p is the chemical po- 
tential, V is the hybridization matrix element (assumed 
k independent), U is the on-site Coulomb interaction, 
and n\ a — / i o ./» )0 - This model describes two electronic 
subsystems, the d and /—electron systems, which are 
coupled by the hybridization matrix. Since the /—states 
form very narrow bands, the Coulomb interaction has a 
strong effect on the /—states. 

It is well known that, at the mean field level, the 
PAM presents a two-band structure, separated by a gap. 
Therefore, the model describes a Kondo-insulator for 
n = 2, describing a metal for other values of n. The 
full density of states shows a sharp peak of width V 2 p^ 
(with po the d— electron density of states) close to the 
renormalized /— level energy, which is responsible for the 
quasiparticle heavy-masses when the chemical potential 
lies within this peak. 

In our study we consider the non-magnetic phase such 
that (nj) = (n^). In the strong coupling limit (U = oo), 
double occupancy of the / sites is forbidden. The sim- 
plest implementation of the U — oo condition is due 
to Coleman,E3 in which the empty /-site is represented 
by a slave boson 6, and the physical operator in 
equation (||) is replaced with 6^/;. Condensation of 
the slave-bosons can be described by the replacement 

h -> {h) = {b\) = y/i. 

The density of the boson condensate z minimizes the 
free energy of the system and the renormalized et en- 
ergy level is obtained after imposing local particle (bo- 
son+fermion) conservation at the /-sites. Using the re- 
sults of Ref. [l3| for the non-superconducting case, the 
mean field equations can be written in terms of the 
Fourier transform of the Green's functions 



G S Ak,T-T') = -{T T h a {T)fl(T')), 
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where N s denotes the number of lattice sites, e/ is the 
renormalized energy of the / orbitals due to the on-site 
repulsion and f(E) is the Fermi function. Equation (||) 
states that the mean number of electrons at an /-site is 
1 — z. For a given number of particles per site, n, these 
equations must be supplemented with the particle con- 
servation condition which yields the chemical potential p 
for any temperature: 



n = 1 — z ■ 
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In the above equations the energies E± and the coherence 
factors u±(e) are given by 

E± = \{e i + e f )±\E{e % ), 



= ~ 2 
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and po(e) represents the density of states of the free 
d— electrons, which we choose of the form 



Po(e) 
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-D < e < D . 
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where D is half the d-electrons bandwidth. 

At weak coupling the on-site Coulomb repulsion term 
Hjj is treated using the usual Hartree-Fock decomposi- 
tion H u — > U (n/) /2 a n { . Since the total number of 
electrons, n^+nf, is fixed, we obtain a mean field equa- 
tion for n f , which is solved together with the equation for 
the chemical potential p. The Coulomb interaction con- 
tributes to the renormalization of eo, leading to a renor- 
malized /— energy given by e/ = eo + U (r?/) /2 . 

Once the mean field equations are solved, the specific 
heat is computed as 



C 



1 d < H > 
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At weak coupling, the mean field parameters are rif and 
p, whereas at U = oo these are rif, e/ and p. In what 
follows we take D — 6 and U, eo, V and T are measured 
in units of D. 

In Figures |l| and |2| we show the specific heat of the 
PAM in the U — oo limit for several values of eo and V, 
respectively. Since we are in the metallic regime the spe- 
cific heat is linear at low temperatures (check inset in Fig. 
|]J). Both the 7 coefficient and the zero-temperature spin 
susceptibility are much higher than the corrresponding 
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values for non- hybridized (i-electrons, signaling the pres- 
ence of heavy quasi-particles. We see that as V or e de- 
crease (with e < 0) the low temperature peak is shifted 
to lower and lower temperatures, whereas the tempera- 
ture position of the high temperature peak remains es- 
sentially constant. Let us define a characteristic temper- 
ature T* as the temperature value where d{d\ a fi,tr)JdT 
has its first local maximum, as can be seen in Fig. [3| for 
two values of V. A similar definition has been adopted 
for the Kondo temperature Tjc (on the lattice) in Ref. 
^. Comparing the results for T* in Fig. ^ with those 
in Figs. |j] and || for C(T), we see that the tempera- 
ture position of the low-temperature peak coincides with 
T* . The dependence of T* on V and eo follows the same 
trends as those found for Tk in Ref. |8|, for the half-filled 
symmetric case, where a Kondo intra-site interaction of 
the form Jk ~ V 2 U/(eo(eo + U)), was assumed to be 
present. Since we have not included the fluctuations of 
the boson fields, it is remarkable that mean field the- 
ory (in the boson fields) can capture this effect since the 
Kondo interactions are not generated at mean field level. 
We interpret the appearence of a specific heat peak at 
temperature T* as being associated with the existence of 
a peak in the density of states. 

As we have seen, in the strong coupling limit U = oo 
the constraint rif < 1 introduces boson fields. At the 
mean field level, the effect of these fields is to renormal- 
ize the hybridization matrix V to V z = \fzV . This im- 
plies that for the same parameters as in the weak cou- 
pling (Hartree-Fock) regime the low temperature peak in 
C(T) shifts to lower temperatures. This can be clearly 
seen comparing Fig. |l| with Fig. pL and Fig. || with Fig. 
||. We have also checked that the low-temperature peaks 
for different V or eo tend to collapse into a single peak 
if T is scaled by T*, as observed for the symmetric case 
in Ref. B, ft is also clear from Figs. |l] and ^| that as 
we move the system towards rif — > 1, which in a more 
elaborated description would correspond to the Kondo 
limit, the single peak structure separates into two peaks. 
It is also clear that the low-temperature peak is in fact a 
superposition of two peaks not completely separated. In 
Ref. H the complete separation of these two peaks took 
place only at very negative values of eo- 

In the single impurity case the model has been solved 
by Bethe Ansatz.E-3 In that case, C(T) also presents a 
double peak structure. In that solution the charge and 
spin excitations are described by different rapidities be- 
coming clear that the low-temperature peak in C(T) is 
associated with the spin excitations whereas the high- 
temperature peak is due to the charge excitations. An 
exact solution is not available for the lattice problem, 
however. 

In more sophisticated treatments of the PAmIIII a lat- 
tice Kondo temperature Tk has been identified with the 
first maximum in the specific heat leading the authors 
to conclude that the first peak is associated with the 
Kondo shielding (spin excitations) whereas the second 



peak is associated with the charge transport (at an en- 
ergy of the order of the hopping taken as unity in this 
paper). But the calculation we just presented shows that 
both the slave-boson mean field theory and the most sim- 
ple Hartree Fock approximation lead to a multiple peak 
structure in the specific heat. These approximations do 
not describe the dynamics or screening of local moments. 
We therefore conclude that the double peak structure in 
the specific heat is simply a band-shape effect, namely, 
the appearence of a peak in the density of states. In- 
deed, such a feature is also present in more sophisticated 
treatments of the PAM.Effl 

If the system undergoes a magnetic or superconduct- 
ing transition, the peak structure is replaced, in general, 
by a lambda-shaped peak signalling the transition to the 
ordered phase or phases if they coexist. In.-the partic- 
ular case of a superconducting ground stateES and zero 
magnetization we have observed that in some cases the 
low-temperature peak is still present. A iStudy of the 
ordered phases will be presented elsewhere.Ej 

In many heavy-fermion materials the /—electrons have 
large total angular momentum. At mean field level and 
for N > 2 the z parameter has a critical temperature 
given by T z = (e/ — /J,)/ ln(iV— f )p,Ja the "Kondo regime" 
6f — fi is a very small quantity.LrE3 As the temperature 
rises V z vanishes. Therefore, we expect for C(T), and 
within mean field theory, a different temperature behav- 
ior from that reported here. 

In summary, we have found that C(T) for the PAM 
presents a double peak structure both in the strong cou- 
pling limit and in the Hartree-Fock approach. This mul- 
tiple peak structure in C(T) is associated with a peak in 
the density of states of the system, due to hibridization 
between delocalized d~ and localized /— states. The 
temperature T* of the first peak corresponds to a maxi- 
mum in the derivative of (d ir7 fi t<7 ) with respect to tem- 
perature. T* decreases as both V and eo decrease. This 
effect is much more pronounced in the strong coupling 
regime, since V is renormalized to smaller values due to 
the factor y/z. 
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FIG. 2. Specific heat as function of T/D for different 
values of V, in the U = oo regime. The parameters are 
eo = -0.57D and n = 1.1. 
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FIG. 1. Specific heat of the PAM in the U = oo limit as 
function of T/D for different values of e . The inset shows a 
linear low-temperature behavior of C(T), as it should be for 
a metal. The parameters are V = 0.17D, n = 1.1. 



FIG. 3. Temperature dependence of the temperature 
derivative of matrix element (dJ ./i i <T for two different values 
of V, in the U — oo regime, for eo = — 0.57D and n = 1.1. 
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FIG. 4. Dependence of the correlation temperature T* as 
a function of: (a) eo and (b) V in the U = oo regime and for 
n = 1.1. 
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FIG. 5. Specific heat of the PAM in the Hartree Fock ap- 
proximation as function of T/D for different values of eo- The 
inset shows a linear low-temperature behavior of C(T), as 
it should be for a metal. The parameters are U = 0.47D, 
V = 0.17D, n = 1.1. 
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FIG. 6. Specific heat of the PAM in the Hartree Fock ap- 
proximation as function of T/D for different values of V. The 
parameters are U = 0.47D, eo = — 0.83D, n = 1.1. 
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